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I.  Introduction 


Early  in  our  attempt  to  create  a  program  for  fitting  the  modified  point-mass  (MPM) 
equation  to  radar  data,  we  noted  a  peculiar  feature  of  this  equation  in  its  present  form: 
it  is  not  given  as  an  explicit  expression  for  the  highest  order  derivative  (in  this  case,  the 

acceleration  vector  U).  Instead,  it  has  the  form 

t  =  f(t,p,u, if)  (i) 

where  P  is  a  parameter  set. 

This  feature  of  U  co-existing  on  both  sides  of  (1)  is  not  an  intolerable  burden  when 
the  aim  is  merely  to  solve  the  equation  numerically,  using  a  specified  set  of  parameter 
values.  Indeed,  for  a  given  set  P,  our  Firing  Tables  Branch  routinely  solves  (1)  by  a 

predictor-corrector  technique  ideally  suited  to  handling  that  ubiquitous  U  (although  addi¬ 
tional  calculations  are  required  at  the  start  of  the  integration). 

In  our  fitting  problem,  however,  the  trajectory  is  known  and  we  want  to  determine  the 
parameter  values  that  yield  the  best  fit  (say,  in  some  least-squares  sense)  to  that  observed 
trajectory.  In  order  to  use  one  of  the  packaged  routines  for  accomplishing  this  curve  fitting, 
we  decided  that  it  was-if  not  absolutely  essential-then  at  least  extremely  desirable  to  work 
with  an  explicit  expression  for  the  acceleration. 

In  the  hope  that  such  a  re-formulation  may  be  of  interest  to  a  wider  audience  than 
the  one  enthralled  by  our  ongoing  trajectory  fitting  problem  (the  latter  group  being  a  very 
small  one  indeed),  we  decided  to  release  the  transformation  as  an  independent  result. 


II.  Modifying  a  Modified  Equation 

The  point  mass  equation  of  motion  can  be  written  in  the  form 

U  =  —AV  +  G  (2) 

where 

U  —  projectile  velocity  with  respect  to  the  earth 
V  =  U  —  W  — -  projectile  velocity  with  respect,  to  the  air 
W  —  wind  velocity  with  respect  to  the  earth 
G  =  the  sum  of  the  gravity  and  Coriolis  accelerations 
and  where 
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(All  symbols  are  defined  in  the  List  of  Symbols.) 

The  MPM  equation  -  the  modified  version  of  Eq.(2)  -  is  derived  in  Reference  l.1  This 
MPM  equation  differs  from  the  point-mass  equation  in  that  the  modified  version  contains 
additional  terms  involving  the  yaw  of  repose,  at.  (The  yaw  of  repose  is  the  steady-state 
angle  of  attack  resulting  from  the  gravity-induced  curvature  of  the  trajectory.)  In  the 
expression  given  in  Ref.  1  for  ae,  we  will  ignore  the  Magnus  moment  coefficient;  that  is, 
we  set 

CMpa  -  0.  (3) 

Under  this  minor  restriction,  the  expression  for  ac  [Eq.(3.3)  of  Ref.  lj  reduces  to 

ae  =  B  U  x  V  (4) 


where 


B  = 


and  where  <f>  is  the  axial  spin  rate. 

Then  the  MPM  equation  [Eq.(3.1)  of  Ref.  1]  can  be  written  in  the  form 

C  =  -AV  +  G  +  M  _  baljLL 


(5) 


where 

is  =  fi'/B  =  V  x  V 

and  where  hi  and  ft  ai  are  dimensionless  coefficients  related  to  the  lift  and  Magnus  forces, 
respectively: 


hL  = 


hM  — 


Cm. 


For  spin-stabilized  artillery  rounds,  MPM  trajectories  based  on  Eq.(5)  -  compared  with 
point-mass  trajectories  -  are  in  much  closer  agreement  with  complete  6DOF  rigid-body 
computer  runs. 

Eq.(5)  differs  from  the  point-mass  equation  (2)  not  only  in  the  addition  of  the  hi  and 
hsi  terms,  but  also  in  the  fact  that  the  drag  coefficient  Cp  contained  in  coefficient  A  is 
now  allowed  to  vary  linearly  with  j <,  | 2 : 

CD  =  CDo  +  CD2  |a,j2.  (6) 

’  “Equations  of  Motion  for  a  Modified  Point  Maas  Trajectory,"  Robert  F.  Lieske  and  Mary  L.  Reiter,  US  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  Report  R1314,  March  1966,  AD  485869. 
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Thus  we  see  that  in  the  MPM  equation  (5),  U  appears  not  only  on  the  left-hand  side  but 
in  every  term  except  G  on  the  right-hand  side.  The  purpose  of  this  report  is  to  present  the 
following  explicit  expression  for  U: 

jj  =  (C  -  A)V  +  gJ]  (7) 

where 


c 


G.4 


ha(G  •  V) 


(1  +  ha)V2 

1 


1  +  ha 


G  + 


hL(G  x  V) 
(1  ~hM)V 


1  —  hm 


—  h\i . 


The  derivation  of  this  form  of  the  MPM  equation  is  given  in  the  Appendix. 

The  attentive  reader  may  be  uneasy  about  that  coefficient  A  in  Eq.(7).  It  does,  after 

all,  contain  Co,  which  in  turn  contains  |ae|.  which  would  seem  to  involve  U.  However, 
— ♦  — * 

a  U- free  expression  for  of  can  now  be  obtained  by  substituting  for  U  from  Eq.(7)  in  the 

definition  of  ac,  Eq.(4).  We  obtain 


n,  r---  BGa  x  f] 

an  expression  devoid  of  V . 


(8) 


III.  Concluding  Remarks 

Note  that  the  appearance  of  U  on  the  right-hand  side  of  the  original  MPM  equation 
(5)  is,  in  a  sense,  illusory.  If  we  use  the  re-defined  yaw  of  repose  from  Eq.(8),  then  (5),  like 

(7),  represents  an  explicit  expression  for  U . 

The  fact  that  the  drag  term  vanishes  in  the  cross-product  U  x  V  is  what  allows  us 

to  declare  Cd  a  function  of  ae  and  still  obtain  an  explicit  expression  for  U.  If  an  unkind 
bailistician  insisted  that  one  or  both  of  the  aerodynamic  parameters  hi  and  must 
depend  on  the  yaw  of  repose  then  no  such  explicit  expression  could  be  obtained. 
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Appendix  :  Derivation  ofEq.(7) 


To  obtain  Eq.(7)  of  the  text,  we  need  a  series  of  preliminary  relationships.  First,  we 
note  that  Eq.(4)  implies  that  ae  -  and  hence  0  -  are  perpendicular  to  V.  Therefore 


0»V  =  0. 


Next  we  note  that  by  Eq.(5), 


~  ^  .  A  .  hL0  hM0xV\  - 

^•V  =  [~AV  +  G  +  ~v - yi  J  *V- 

But  the  hi  term  above  is  zero  by  Eq.(Al)  and  the  term  is  zero  by  the  triple  scalar 
product  rule.  Hence 

V  •  V  =  -AV2  +  G»V  =  hb  (A2) 

where  the  symbol  hb  is  introduced  as  a  temporary  convenience.  Then 

0  x  V  =  (U  x  V)  x  V 

=  hbV  -  V2U  (.43) 

by  the  triple  vector  product  rule.  This  same  rule  gives  us 

{0  xV)x  V  =  (0  •  V)V  -  V20 

=  -V20.  (.44) 

We  can  now  obtain  a  convenient  expression  for  0: 


0  =  V  x  V 


—  AV  4  G  4- 


hL0  hM0  x  Vr 


,-r  r.  hL{0xV)  hM(0xV)xV 
=  6xV  +  — - - V3 - 


so  that  by  Eqs.(A3)  arid  (A4), 


3  .  >>dhV-\nV)  '  ,  s 

p  =  G  x  \  + - - - +  xMii 


Combining  the  0  terms,  we  obtain 

s  Gy  v +  (hMV)V  ~(hLV)U 
1  -  hM 

Substituting  (A3)  and  (A5)  in  Eq.(5)  of  the  text,  we  obtain 


on 


v  -  (AUOt'l  -  [m7  -  V'O 
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Combining  similar  terms,  we  have 


1  + 


f‘l 


1  -  h 


At 


—  !>  m  I  t  = 


-A  + 


ht.  (  h\ 


l'}  1  -h 


- ^Af 

M 


V  +  G  + 


or 


But  from  Eq.(A2). 


,  ft  Jih  (  G  •  V 

~A+-rrr  =  ~A  +  ha  -.4  + 


r2 


—  -  ( 1  +  ha  )A  + 


V2 

haG  •  V 
V2 


Substituting  (AT)  in  (AG),  we  have 


h^Lf  X  1 


(A6) 


(AT) 


( 1  +  hQ  )C  =  -( 1  +  K  )AY  +  (  -'4^  I  V  +  C  +  /l,  c  X  r 


r2 


(1  -  /«Af)V' 


from  which  Eq.(T)  of  the  text  follows  at  once. 
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List  of  Symbols 


A  (SCo)  (V) 

r  .hs&S-l 

(1  +ha)V* 

Cp  drag  coefficient,  |dra<7  /orcej/(pV25/2) 

Cdo ,  Cp?  zero-yaw  and  yaw-drag  coefficients: 

Cp  —  Cp0  +  Cp2  |oe|2 

C(,  roil  damping  moment  coefficient: 

|ro//  damping  momcnt\  =  ±(pV2 S£ /2)(4>(! /V)C/p 

CLa  lift  force  coefficient: 

| lift  forcc\  =  ±(pV2 S /2)\ae\C L.a 

CAf  Magnus  moment  coefficient: 

\Magnus  moment]  —  ±(pV2SC/2)(4>C/V)\ot\CMf>a 

C\ia  static  moment  coefficient: 

js/a/fc  moment j  =  ±(pV25t’/2)|cT(.|CA/0 


Cv, 

G 

Ga 

ha 

hb 

h, 

h\i 


Magnus  force  coefficient: 

j  Magnus  force  |  -  ±(pV2  S/2)(<j>C/V)\ae\CNp 

the  sum  of  gravity  and  Coriolis  accelerations 
u  i  (1-AW)V 


_ 1  _ 

l+'.a 


l-h* 


-  h 


M 


U»Y  =  -AV2  +  G»V 

«(§fc)(f) 

y (tr)  t  *0’ 

axial  moment  of  inertia 
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t 

m 

MPM 

R 

S 

t 

u 


V 

V 
w 

Qe 

P 

p 

4> 

(') 


Ix  /  mC2 

reference  length 
projectile  mass 
Modified  Point  Mass 

effective  radius  of  the  earth  (6  356  766  m) 

reference  area,  7r£2/4 

time 

projectile  velocity  with  respect  to  the  earth 

m 

U  —  W,  projectile  velocity  with  respect  to  the  air 
wind  velocity  with  respect  to  the  earth 

the  yaw  of  repose,  B  U  x  V  =  B  GA  x  V 

S'/B 

air  density 
axial  spin  rate 
d(  )/dt 
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